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^0 I Gravitomagnetic field is covariantly split in the intrinsic and extrinsic parts, which 



are generated by rotational and translational currents of matter respectively. The 



Qh< intrinsic component has been recently discovered in the LAGEOS spacecraft experi- 

\ ment. We discuss the method of detection of the extrinsic tidal component with the 

lunar laser ranging (LLR) technique. Analysis of the gauge residual freedom in the 



relativistic theory of three-body problem demonstrates that LLR is currently not 
capable to detect the extrinsic gravitomagnetic effects which are at the ranging level 
of few millimeters. Its detection requires further advances in the LLR technique that 
are coming in the next 5-10 years. 



> 
(N 

0\ . PACS numbers: 04.20.-q, 04.80. Cc, 96.25.De 

m 
cn 

^ ■ Detection of gravitomagnetic components as predicted by Einstein's theory of general rel- 
^ \ ativity is one of the primary goals of experimental gravitational physics. The paper |l| 
states that the gravitomagnetic interaction plays a part in shaping the lunar orbit readily 
obervable by LLR. The authors picked up a "gravitomagnetic" term from the parameterized 
post-Newtonian (PPN) equation of motion of massive bodies [2;] and proved that it correctly 



reproduces the Lense-Thirring precession of the GP-B gyroscope. The paper |l| argues that 
the very same term in the equations of motion of the Moon, derived in the solar-system 
barycentric (SSB) frame, perturbs the lunar orbit with a radial amplitude ~ 6 meters that 
was observed with the lunar laser ranging (LLR). We have explained and confirm in more 
detail in the present paper that the gauge freedom of the relativistic three-body problem 
suppresses the gravitomagnetic effects in the lunar motion, that depend on the Earth's ve- 
locity V around the Sun, to the level < 1 centimeter. It makes LLR currently insensitive to 
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the gravitomagnetic interaction. 

There are two types of mass currents in gravity 0, [sl- The first type is produced by the 
intrinsic rotation of matter around body's center of mass. It generates an intrinsic gravito- 
magnetic field tightly associated with body's angular momentum (spin) and most research 
in gravitomagnetism has been focused on the discussion of its various properties. Textbook 
6| gives a comprehensive review of various aspects of the intrinsic gravitomagnetism. It 
is interesting to note that the intrinsic gravitomagnetic field can be associated with the 
holonomy invariance group [tI . Some authors [s], [qJ have proposed to measure the intrin- 
sic gravitomagnetic field by observing quantum effects of coupling of fermion's spin with 
the angular momentum of the Earth. It might be worthwhile to explore association of the 
intrinsic gravitomagnetism with the classic Hannay precession phase 1128. 

The first classic experiment to test the intrinsic gravitomagnetic effect of the rotating 
Earth has been carried out by observing LAGEOS in combination with other geodetic satel- 
lites , ll5|, llGj which verified its existence with a remarkable precision as predicted by 
Einstein's general relativity. Independent experimental measurement of the intrinsic gravit- 
omagnetic field of the rotating Earth is currently under way by the Gravity Probe B mission 



17|. 



The second type of the mass current is caused by translational motion of matter. It 
generates an extrinsic gravitomagnetic field that depends on the frame of reference of ob- 
server and can be either completely eliminated in the rest frame of the matter or significantly 
suppressed by the transformation to the local inertial reference frame of observer. This prop- 
erty of the extrinsic gravitomagnetic field is a direct consequence of the gauge invariance of 
Einstein's gravity field equations for an isolated astronomical system 18| embedded to the 
asymptotically-flat space-time. Experimental testing of the extrinsic gravitomagnetic field 
is as important as that of the intrinsic gravitomagnetic field. The point is that both the 
intrinsic and the extrinsic gravitomagnetic fields obey the same equations and, therefore, 
their measurement would essentially complement each other [s]. Furthermore, detection of 
the extrinsic gravitomagnetic field probes the time-dependent behavior of the gravitational 
field which is determined by the structure of the gravity null cone (the domain of causal 
influence) on which the gravity force propagates. Experimental verification of the gravito- 
magnetic properties of gr avity is important for the theory of brarrewor.ds Q arrd for setting 
other, more stringent limitations on vector-tensor theories of gravity [20] . 
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Ciufolini ^ proposed to distinguish the rotationally-induced gravitomagnetic field from 
the translationally-induced gravitomagnetic effects by making use of two scalar invariants 
of the curvature tensor 



/2 = R^pf^^R^f^P^'E^ 



fiupa } 



(1) 

(2) 



where Ra/Sfiu is the curvature tensor, E^yp^ is the fully anti-symmetric Levi-Civita tensor with 



det((7^,y) < is the determinant of the metric tensor. Ciufolini 



^0123 = +v^' and g 

notices that a weak gravitational field of an isolated astronomical system yields /2 = if the 
intrinsic gravitomagnetic field is absent. However, one should not confuse the invariant I2 
with the gravitomagnetic field itself. The gravitomagnetic field is generated by any current 
of matter. Hence, J2 = does not mean that any gravitomagnetic field is absent as has 



been erroneously interpreted in 



2l|. Equality /2 = only implies that the gravitomagnetic 



field is of the extrinsic origin (Ji 7^ 0), that is generated by a translational motion of 
matter. The translational gravitomagnetic field can be measured, for instance, by observing 
the gravitational defiection of light by a moving massive body like Jupiter 22l, 1231]. This 
gravitomagnetic frame-dragging effect on the light ray was indeed observed in a dedicated 
radio- interferometric experiment 241 1. 

n 

Paper [1[ makes an attempt to demonstrate that the extrinsic gravitomagnetic field can 
be measured by making use of the LLR observations of the lunar orbit. This must not be 
confused with the measurement of the intrinsic gravitomagentic field by means of the satellite 
laser ranging technique applied to LAGEOS Q, Q, Q Q. The LAGEOS experiment 
measures the Lense-Thirring precession of the satellite's orbit caused by the Earth's angular 
momentum entering g^i component of the metric tensor. The authors of [l| have been trying 
to measure the gravitomagnetic precession of the lunar orbit caused by the orbital motion 
of the Earth-Moon system around the Sun. They used the barycentric coordinates of the 
solar system (BCRS) to derive the equation of motion of the Moon relative to the Earth. 
The equation is effectively obtained as a difference between the Einstein-Infeld-Hoffmann 
(EIH) equations of motion for the Earth and for the Moon with respect to the barycenter 
of the solar system, and repeats the original derivation by Brumberg 25|, l26|] , which later 



was independently derived by Baierlein 



27|. 



The barycentric equation of motion of the Moon formally includes the gravitomagnetic 



perturbation in the following form 



7-1 

o-GM = a H , (3) 



where 7 parameterizes a deviation from general relativity, the bold letters denote spatial 
vectors, and the dot between two vectors means their Euclidean dot product. The general- 
relativistic post-Newtonian acceleration 

a=^[r{V-u)-V{u.r)] , (4) 

where m is mass of the Earth, r is radius of the lunar orbit, r is the unit vector from the 
Earth to the Moon, V is the Earth's velocity around the Sun, and u is the Moon's velocity 
around the Earth. 

We notice ^ that the barycentric coordinate frame referred to the geocenter by a simple, 
Newtonian-like spatial translation (the time coordinate is unchanged) 

r = X- XE{t) , (5) 

as it is obtained in 



is not in a free fall about the Sun, and does not make a local inertial 
frame. Thus, perturbations in Eqs. ©-(jlj) can not be interpreted as physically observable 
and, in fact, represent a spurious gauge-dependent effect that is canceled by transformation 
to the local-inertial frame of the geocenter. This transformation is a generalized Lorentz 
boost with taking into account a number of additional terms due to the presence of the 
external gravitational field of the Sun 13 

The gauge freedom of the lunar equations of motion must be analyzed to eliminate all 
gauge-dependent, non-observable terms. Only the terms in the equations of motion, which 
can not be eliminated by the gauge transformation to the local inertial frame can be phys- 
ically interpreted. The analysis of the gauge freedom in the three body-problem had been 
done in 



30 



31 



32l |. It proves that all non-tidal and V-dependent terms, including the first 
term in the right side of Eq. ([3]), are pure coordinate effects that disappear from the lunar 
equations of motion after transformation to the geocentric, locally-inertial frame. This is 
because the Lorentz invariance and the principle of equivalence reduce the relativistic equa- 
tion of motion of the Moon to the covariant equation of the geodesic deviation between 



the Moon's and the Earth's world lines 



31 



32l |. where gravitomagnetic effects appear only 



as tidal relativistic forces with amplitude smaller than 1 centimeter. The covariant nature 



of gravity tells us that if some effect is not present in the local frame of observer, it 
can not be observed in any other coordinate system. This means that besides physically- 
observable terms, the barycentric LLR model jl|, 13] also operates with terms having the 
gauge-dependent origin, which mathematically nullify each other in the data-processing com- 
puter code irrespectively of the frame of reference. The mutually annihilating terms enter 
different parts of the barycentric LLR model with opposite signs 3, 3li but, if taken sepa- 
rately, can be erroneously interpreted as really observable. This is what exactly happened 
with the misleading analysis given in 

General relativity indicates that the barycentric EIH lunar equations of motion may admit 
the observable gravitomagnetic acceleration only in the form of the second term in the right 
side of Eq. ((31) that is proportional to 7 — 1. Radio experiments set a limit on 7 — 1 < 10~^ ^ 
that yields IcigmI < 1 millimeter. The current half-centimeter accuracy of LLR is insufficient 
to measure such negligible effect. Moreover, our analysis of the gauge invariance of the scalar- 



tensor theory of gravity [30|] points out that the term being proportional to 7 in equation 
([3]) is also eliminated in the locally-inertial, geocentric reference frame. We conclude that 
LLR is currently insensitive to the gravitomagnetism and, yet, can not compete with the 
LAGEOS and/or GP-B experiments. 



Recent paper by Soffel et al. [33| is another attempt to prove that the extrinsic grav- 
itomagnetic acceleration (jl]) can be measured with the LLR technique in a locally-inertial 
reference frame. The authors of 33| accept our criticism 3| but continue to believe that 



one can measure the extrinsic gravitomagnetic field by making use of the preferred frame 



parametrization of the gravimagnetic terms. To this end, Soffel et al. 33|] introduce the 
preferred-frame generalization of equation ([3]) by replacing 

7-1^7-1 + r/G/4 , (6) 

where rjc is a parameter labeling the gravitomagnetic effects in the barycentric equations of 
motion of the Moon. Parameter rjc = —ai/2 in the framework of the PPN formalism [2I 
but unless this is not stated explicitly one can regard rjc as an independent parameter and 
fit it to LLR data irrespectively of ai. This was done by Soffel et al. [33] who had obtained 
riG = (0.9 ±0.7) X 10~3_ 

We argue that this measurement says nothing about the extrinsic gravitomagnetic field. 
This is because rjc has no any fundamental significance. Its value is not invariant and 
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crucially depends on the choice of the preferred reference frame, which one uses for processing 
LLR data. In general theory of relativity r^^ = in any frame, that means a vanishing (non- 
observable) effect. Indeed, analytic calculations given in appendix of the paper {sH] and in 



34 



35| confirm that in the PPN framework with ai = all gravitomagnetic effects in the 



motion of the solar system are nullified by the corresponding gravitoelectric effects from 
other well-established post-Newtonian gravitational potentials. Those calculations reveal 
that the gravitomagnetic effect described by equation ([3]) is nothing but a symbolic property 
of the particular coordinate system used for calculations. LLR data fit and the coordinate- 
dependent limits on t]g obtained in paper 33|] confirm that the gravitational model used 
in the data processing, is self-consistent. However, it neither means that the extrinsic 
gravitomagnetic field was measured nor that general theory of relativity was tested. This 
is because the gauge-invariance is the main property of a large class of the metric-based 
gravitational theories and testing the self-consistency of the LLR equations does not single 
out a specific gravitational theory. 

In order to measure the extrinsic gravitomagnetic field, one has to find the real grav- 
itomagnetic effect in the motion of the Moon, which does not vanish in the framework of 
general theory of relativity in the locally- inert ial frame of observer on the Earth. For this 
reason, the primary goal of the relativistic theory of the lunar motion is to construct an 
inertial reference frame along the world-line of the geocenter and to identify the gravito- 
magnetic effects in this frame. This task was solved in our paper [31] in the post-Newtonian 
approximation in the case of a three-body problem (Sun, Earth, Moon) under assumption 
that the Moon is considered as a test particle. In the quadrupole approximation the metric 
tensor in the locally-inertial geocentric reference frame X" = {cT, X) has the following form 



Goo(T,X) = -1 + - 



U{T,X)+QpX^ + ^Qp,XPX'' 



+ 0{-,], (7) 



Go^iT,X) = [W{T,X) + eipkCp,X''X'^]+o(^^^ 



2 

Gij(T, X) = 6ij H — - 



c2 



UiT,X)+Q,XP + ^Q,gXPX'^ 



(8) 

+ ^ (^) ' 



where U{T,X) is the Newtonian potential of the Earth, W{T,X) is the post-Newtonian 
vector-potential of the Earth, Qp is the acceleration of the geocenter with respect to the 
geodesic world line, Qpg and Cpg are tidal gravitational-force gradients from the Sun, and 
eipk is the fully anti-symmetric symbol. 
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We notice that W{T,X) represents the intrinsic gravitomagnetic field of the Earth, 



36| that 



which has been measured in the LAGEOS experiment [l^. It can be shown 
gravitomagnetic vector potential U^{T,X) produces a negligibly small acceleration on the 
lunar orbit, and can be discarded. The tensor potential Cpg has the extrinsic gravitomagnetic 
origin as it is generated in the geocentric frame by the motion of an external mass (the Sun) 
with respect to the Earth. This potential is expressed in terms of the orbital velocity of the 



Earth, V^, and the Newtonian tidal matrix Q 



pq 



Cpg = eikp (^iQkq — VkQiq + -^kqV'-'Qij — -^iqV^Qkj^ , (lO) 

where 

Qpq = {^XpXq — 6pgX^) , (11) 

M is mass of the Sun, and R is the distance between the Sun and Earth. 

The gravitomagnetic tidal force causes a non-vanishing gravitomagnetic acceleration of 
the Moon, which reads in the locally- inert ial geocentric frame as follows (see equation 7.12 
from [31.] ) 

^GM ~ ^'^^ijkCjgU'^r'' , (12) 

where is the geocentric velocity of the Moon, and r* is the Earth-Moon radius-vector. 
This gravitomagnetic acceleration causes a radial oscillations of the lunar orbit that can be 
estimated by making use of equations (8) and (9). Noticing that the term Qijr^ is propor- 
tional to the Newtonian tidal force from the Sun, which produces the variation inequality 
of the lunar orbit [s^l, one obtains 

Vu 

\Acm\ — 10 X (variation) x — ^ . (13) 



The variation amounts to ~ 3000 km in radial oscillation [37]. Hence, the amplitude of the 
oscillation of the lunar orbit caused by the gravitomagnetic tidal force has an amplitude of 
about 1 cm. This is unmeasurable with the current LLR data but can be measured in the 
next 5-10 years with the advent of a millimeter-range LLR technology 38 1. 
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